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A precise definition for a quantum electron thermometer is given, as an electron reservoir coupled 
locally (e.g., by tunneling) to a sample, and brought into electrical and thermal equilibrium with 
it. A realistic model of a scanning thermal microscope with atomic resolution is then developed, 
where the resolution is limited in ultrahigh vacuum by thermal coupling to the electromagnetic 
environment. We show that the temperatures of individual atomic orbitals or bonds in a conjugated 
molecule with a temperature gradient across it exhibit quantum oscillations, whose origin can be 
traced to a realization of Maxwell's demon at the single-molecule level. These oscillations may be 
understood in terms of the rules of covalence describing bonding in 7r-electron systems. Fourier's 
law of heat conduction is recovered as the resolution of the temperature probe is reduced, indicating 
that the macroscopic law emerges as a consequence of coarse graining. 



Recent advances in thermal microscopy^— have opened 
the door to understanding noncquilibrium thermody- 
namics at the nanoscale. The nonequilibrium tempera- 
ture distribution in a quantum system subject to a ther- 
mal or electric gradient can now be probed experimen- 
tally, and a number of fundamental questions can be 
addressed: Can significant temperature variations occur 
across individual atoms or molecules without violating 
the uncertainty principle? How are the electronic and 
lattice temperatures related in a nanostructure out of 
thermal equilibrium? How does the classical Fourier law 
of heat conduction emerge^ from this quantum behavior 
in the macroscopic limit? 

In order to address these questions theoretically, a def- 
inition of a nanoscale thermometer that is both realis- 
tic and mathematically rigorous is needed. According 
to the principles of thermodynamics, a thermometer is 
a small system (probe) with some readily identifiable 
temperature-dependent property that can be brought 
into thermal equilibrium with the system of interest 
(sample). Once thermal equilibrium is established, the 
net heat current between probe and sample vanishes^ - — 
and the resulting temperature of the probe constitutes 
a measurement of the sample temperature. High spatial 
and thermal resolution require that the thermal coupling 
between probe and sample is local^ and that the cou- 
pling between the probe and the ambient environment is 
smallji respectively 

It should be emphasized that, out of equilibrium, the 
temperature distributions of different microscopic de- 
grees of freedom (e.g., electrons and phonons) do not, in 
general, coincide, so that one has to distinguish between 
measurements of the electron temperature^^ and the 
lattice temperature.— ^ This distinction is particularly 
acute in the extreme limit of elastic quantum transport,-^ 
where electron and phonon temperatures are completely 
decoupled. 

In this article, we develop a realistic model of a scan- 
ning thermal microscope (SThM) operating in the tun- 



neling regime in ultrahigh vacuum, where the vacuum 
tunneling gap ensures that phonon heat conduction to 
the probe is negligible. Since electrons carry both charge 
and heat, an additional condition is necessary to define 
an electron thermometer. We proceed by noting that, as 
a practical matter, in order to reduce the thermal cou- 
pling of the thermometer to the ambient environment, it 
should form an open electrical circuit (or have very high 
impedance to ground). This ensures that, in addition 
to the heat current, the electrical current between sam- 
ple and probe vanishes. This is also in accord with the 
common-sense notion that thermometers are not current 
sources or sinks. An electron thermometer is thus de- 
fined as an electron reservoir whose temperature is fixed 
by the conditions of electric and thermal equilibrium with 
the sample: 

7^=0,^ = 0,1, (1) 

where —elp ^ and ij^ are the electric current and heat 
current, respectively, flowing into the probe. In an ideal 
measurement, the sample would be the sole source of 
charge and heat flowing into the probe, but we also con- 
sider nonideal measurements, where there is an addi- 
tional thermal coupling to the ambient environment. In 
practice, this coupling plays a crucial role in limiting the 
resolution of temperature measurements i 

In a measurement of the temperature distribution in a 
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FIG. 1: Schematic representation of a temperature probe as 
the third terminal of a thermoelectric circuit. 
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conductor subject to thermal and/or electric gradients, 
the electron thermometer thus serves as the third termi- 
nal in a three-terminal thermoelectric circuit, a gener- 
alization of Buttiker's voltage probe concept^ 3 , (see Fig. 
[1]). Note that the conditions ([1]) allow a local temper- 
ature to be defined under general thermoelectric bias 
conditions, relevant for the analysis of nonequilibrium 
thermoelectric device performance^^— Previous theo- 
retical analyses of quantum electron thermometers cither 
completely neglected thermoelectric effects^ or consid- 
ered the measurement scenario ([1} as only one of several 
possibilities ; 18 ' 19 while the subtle definition of local tem- 
perature given in Ref. [l^ (thermometer causes minimal 
perturbation of system dynamics) may capture the spirit 
of our two separate conditions, at some level. A recent 
review of the topic is given in Ref. Ho- 
using our model of a nanoscale electron thermometer, 
we investigate the nonequilibrium temperature distribu- 
tions in single-molecule junctions subject to a thermal 
gradient. Quantum temperature oscillations analogous 
to those predicted in one-dimensional systems^ are pre- 
dicted in molecular junctions for several different con- 
jugated organic molecules, and are explained in terms of 
the rules of covalence describing bonding in 7r-conjugatcd 
systems. In terms of directing the flow of heat, the rules 
of covalence can be seen as an embodiment of Maxwell 's 
Demon at the single-molecule level. 

It has been argued that in some systems, quan- 
tum temperature oscillations can be washed out by ei- 
ther dephasing£ or disorder^ leading to restoration of 
Fourier's classical law of heat conduction. However, in 
molecular junctions the required scattering would be so 
strong as to dissociate the molecule. We investigate 
the effect of finite spatial resolution on the nonequilib- 
rium temperature distribution, and find that Fourier's 
law emerges naturally as a consequence of coarse-graining 
of the measured temperature distribution. Thus our res- 
olution of the apparent contradiction between Fourier's 
macroscopic law of heat conduction and the predicted 
non-monotonic temperature variations at the nanoscale 
is that the quantum temperature oscillations are really 
there, provided the temperature measurement is carried 
out with sufficient resolution to observe them, but that 
Fourier's law emerges naturally when the resolution of 
the thermometer is reduced. 



I. ELECTRONIC TEMPERATURE PROBE 

Consider a general system with M electrical contacts. 
Each contact a is connected to a reservoir at temperature 
T a and electrochemical potential fi a . In linear response, 
the electrical current 
into reservoir a may be expressed as 



ela^ and heat current Ia^ flowing 



r(V) 



M 

£ 

13=1 



C 



(H-l) 

a/3 



T 



-(T/3 — T a ) 



(2) 



where is a linear- response coefficient. Eq. @ is a 
completely general linear-response formula, and applies 
to macroscopic systems, mesoscopic systems, nanostruc- 
tures, etc., including electrons, phonons, and all other 
degrees of freedom, with arbitrary interactions between 
them. For a discussion of this general linear-response 
formula applied to bulk systems, see Ref. l2ll 

In this article, we consider systems driven out of equi- 
librium by a temperature gradient between reservoirs 1 
and 2. Thermoelectric effects are included, so the chem- 
ical potentials of the various reservoirs may differ. We 
consider pure thermal circuits (i.e., open electrical cir- 
cuits), for which J a °' ) = OVa. These conditions may be 
used to eliminate the chemical potentials /i a from Eq. 
([2]), leading to a simpler formula for the heat currents 
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In the absence of an external magnetic field C^l = £ 
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and the three-terminal thermal conductances are given 
by 
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This paper is organized as follows: In Sec. U a gen- 
eral linear-response formula for an electron thermometer 
is derived. A realistic model of a scanning thermal mi- 
croscope with sub-nanometer resolution is developed in 
SecHH including a discussion of radiative coupling of the 
probe to the environment. Results for the nonequilib- 
rium temperature distributions in single-molecule junc- 
tions subject to a thermal gradient are presented in Sec. 
IIIII A discussion of our conclusions is given in Sec. IIVI 
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An equivalent circuit for and is given in Fig. [2] 
The first line of Eq. @ resembles the famil- 
iar two-terminal thermal conductanc e) 16 ' 22 ' 23 K Q( g = 
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with £^°2 replaced by Eq. (|5|). 
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FIG. 2: Equivalent circuit for and where the re- 

sistance R a p = [e 2 £^] -1 . e 2 C^ is the loop conductance of 

the circuit and e 2 £^°j is the effective two-terminal conduc- 
tance between terminals a and /3. 

(2) 

Since C a L is usually the dominant term, R a p is often 
comparable to the two-terminal form n a p (cf. Fig. [5]). 
However, the discrepancy is sizable in some cases. Al- 
though it might be tempting to interpret the second line 
in Eq. Q as a nonlocal quantum correction to the ther- 
mal conductance, it should be emphasized that this is a 
generic three-terminal thermoelectric effect, that arises 
in bulk systems as well as nanostructures. 



II. QUANTUM ELECTRON THERMOMETER 

We consider nanoscale junctions with weak electron- 
phonon coupling operating near room temperature. Un- 
der linear-response conditions, electron-phonon interac- 
tions and inelastic scattering are weak in such systems, 
and the indirect phonon contributions to C ®p and C^p 
can be neglected, while the direct phonon contribution 

(2) 

to C p p is zero due to the vacuum tunneling gap. The 
linear response coefficients needed to evaluate Eq. © 
may thus be calculated using clastic electron transport 
theoryl 6 * 2 ^ 

41 = lJ dE ( E - Mof T aP (E) (-§§) , (9) 

where fo is the equilibrium Fermi-Dirac distribution of 
the electrodes at chemical potential and temperature 
To. The transmission function may be expressed a s 12 ' 24 

T a p(E) = Tr {T a (E)G(E)T^E)G\E)} , (10) 

where r° (E) is the tunneling- width matrix for lead a and 
G(E) is the retarded Green's function of the junction. 



A. Temperature Measurement 

In addition to the coupling of the temperature probe 
to the system of interest, we assume the probe also has 
a small thermal coupling k p q to the environment at tem- 
perature To. The environment could be, for example, the 
black-body radiation or gaseous atmosphere surrounding 
the circuit. The heat current flowing from the environ- 
ment into the probe must be added to Eq. §S§ to deter- 
mine the total heat current: 

2 

^^/W^-^+MTo-Tp). (7) 

/3=1 

Thermal coupling to the environment is important when 
the coupling to the system is weak, and is a limiting fac- 
tor in the thermal resolution of any temperature probe. 
The environment is effectively a fourth terminal in the 
thermoelectric circuit, but since its electrical coupling to 
the system and probe is zero, the thermal conductances 
lipi , k P 2 in Eq. ([7]) have the three-terminal form (01 . Solv- 
ing Eqs. (JT|) and (J7JI for the temperature, we find 

kpiTi + k p2 T2 + k p qT 

1 p - — - — — - — ; \°) 

K p l + K p 2 + K p o 

for a probe in thermal and electrical equilibrium with, 
and coupled locally to, a system of interest. Equations 
(U) and © provide a general definition of an electron 
thermometer coupled to a system with a temperature 
gradient across it, in the linear response regime. 



A. SPM-based temperature probe of a 
single-molecule junction 

As an electron thermometer with atomic-scale reso- 
lution, we propose using a scanning probe microscope 
(SPM) with a conducting tip mounted on an insulating 
piezo actuator designed to minimize the thermal cou- 
pling to the environment. The tip could serve e.g. as a 
bolometer or thermocouple^ and its temperature could 
be read out electrically using ultrafine shielded wiring. 
The proposed setup is essentially a nanoscale version of 
the commercially available SThM, and is analogous to the 
ground-breaking SThM with lOnm resolution developed 
by Kim et al4 

Such an atomic-resolution electron thermometer could 
be used to probe the local temperature distribution in 
a variety of nanostructures/mesoscopic systems out of 
equilibrium. In the following, we focus on the specific 
example of a single-molecule junction (SMJ) subject to 
a temperature gradient, with no electrical current flow- 
ing. In particular, we consider junctions containing con- 
jugated organic molecules, the relevant electronic states 
of which are determined by the 7r-orbitals. We consider 
transition metal tips, where tunneling is dominated by 
the d-like orbitals of the apex atom^ 

The tunnel coupling between the tip of the electronic 
temperature probe and the 7r-systcm of the molecule is 
described by the tunneling-width matri x 25 ' 26 

r% m (E) = 2irV n V£ P p(E), (ff) 

where n and m label 7r-orbitals of the molecule, p P {E) 
is the local density of states on the apex atom of the 
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FIG. 3: The calculated spatial map of Tr{r p } for a Pt 
electron thermometer scanned 3.5A above the plane of 
a benzene molecule and a schematic representation of a 
para-benzenedithiol ([1,4]BDT) junction. Peak values of 
Tr{r p } ~16.6meV correspond to the centers of the carbon 
atoms, which are numbered in black. The sulfur and gold 
atoms were drawn using their covalent radii of 102pm and 
134pm, respectively. 
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FIG. 4: Spectral function A{E) = -l/7rimTr{G} of a 
[1,4]BDT (para) junction in the vicinity of the HOMO and 
LUMO resonances, calculated using many-body theory in- 
cluding the effect of the the partially charged sulfur atoms 
on the intramolecular potential. In all benzene simulations 
Ti = F2 =0.69eV and the excess charge on the thiol end- 
groups ~-0.29e, values which give the best agreement with the 
measured thermopower— and linear-response conductance^ 
measurements. The blue line indicates the Fermi energy of 
gold EP=-5.53eV, averaged over the [110], [111], [320], and 
[331] crystal planes.— 



diativc, so that 



probe, and V m is the tunneling matrix element between 
the evanescent tip wavcfunction and orbital m of the 
molecule. Since the temperature probe is in the tun- 
neling regime, and not the contact regime, the phonon 
contribution to the transport vanishes; heat is exchanged 
between system and probe only via the electron tunneling 
characterized by T p . 

In Fig. [31 the trace of T p (Ep) is shown for a Pt 
temperature probe held 3.5A above the plane of a Au- 
[l,4]benzencdithiol-Au ([1,4]BDT) molecular junction. A 
schematic of the [1,4]BDT junction is also shown, with 
sulfur and gold atoms drawn to scale using their covalent 
radii of 102pm and 134pm, respectively. Peaks in Tr {r p } 
correspond to the locations of carbon 7r-orbitals, labeled 
with black numbers. 

The density of states (DOS) of the para BDT junction 
is shown in Fig. 21 simulated using our many-body the- 
ory including the electrostatic influence of the thiol end 
groups. The blue vertical line is set at the Pt Fermi en- 
ergy Ep=-5.53cV averaged over the [110], [111], [320], 
and [331] crystal planes^ Our many-body theory accu- 
rately reproduces the fundamental gap of gas-phase ben- 
zene (~10.4eV) allowing us to unambiguously determine 
the energy-level alignment between the electrodes and 
molecule. Transport occurs within the HOMO-LUMO 
gap, but is dominated by the HOMO resonance. This is 
true for all of the molecular junctions considered in this 
article. 



B. Radiative coupling to the environment 

Finally, we assume the coupling of the temperature 
probe to the ambient environment is predominantly ra- 



n pa = 4eAaTl (12) 

where e and A are the emissivity and surface area, re- 
spectively, of the metal tip, and a = (TT 2 kg/6Qh 3 c 2 ) is 
the Stefan-Boltzmann constant. The tip coupling to the 
environment should be weak compared to the quantum 
of thermal conductance ko = (tt 2 /3)(k 2 3 T/h) in order to 
resolve quantum effects on the temperature within the 
SMJ. For a "spherical tip" with A = 4ttR 2 , 

K p o _ 8tt 2 £ / Rk B T \ 2 Tq=30ok Q27c f R \ / 13 n 
kq 5 \ he J \A™ / 

The conducting tip of the temperature probe must thus 
have linear dimensions R < 1/im in order to resolve quan- 
tum effects at room temperature. A conducting tip of 
small volume will also ensure rapid equilibration of the 
probe with the sample. 

We do not include the direct radiative contribution to 
K p i and k P 2- Since the separation between electrodes 1 
and 2 is much less than the photon thermal wavelength, 
we consider that black-body radiation from the two elec- 
trodes contributes to a common ambient environment at 
temperature To. 

In all simulations presented here, we consider a Pt 
electron thermometer with an effective blackbody sur- 
face area equal to a 50nm sphere with an emissivity 
of 0.1j2£ so that k p o/kq = 6.75 x 10~ 5 . Temperature 
probes constructed from other metals are predicted to 
exhibit qualitatively similar effects. Eq. (|12p may rep- 
resent an overestimate due to quantum suppression of 
radiative heat transfer for structures smaller than the 
thermal wavelength^ in which case the conditions on 
the probe dimensions would be less restrictive. 
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III. RESULTS 

In this section, we investigate the spatial tempera- 
ture profiles for three Au-benzenedithiol-Au (BDT) junc- 
tion geometries, a linear [I,6]hcxatriencdithiol junction 
and a poly-cyclic [2,7]pyrcncdithiol junction. The BDT 
and hexatrienc SMJ calculations were performed using 
a molecular Dyson equation (MDE) many-body trans- 
port theory^ in which the molecular 7r-system is solved 
exactly, including all charge and excited states, and the 
lead-molecule tunneling is treated to infinite order. The 
transport calculations for the pyrene junction were per- 
formed using Hiickcl theory. In all cases, the ambient 
temperature is taken as T Q ~ 300K. 



A. Benzenedithiol junctions 

We investigate the temperature distributions for three 
Au-benzenedithiol-Au (BDT) junction geometries: the 
'para' [1,4]BDT junction, shown schematically along with 
the trace of the lead-molecule coupling matrix in Fig. [3l 
the 'ortho' [1,2]BDT junction, and the 'mcta' [1,3]BDT 
junction. The calculated spatial temperature distri- 
bution for a para junction is shown in Fig. [51 with 
T t = 325A" and T 2 = 275A". The figure illustrates 
quantum oscillations of the local temperature near the 
molecule, which are clearly resolvable using our model of 
a nanoscale electron thermometer. Each of the 7r-orbitals 
of the molecule has a characteristic temperature diffcr- 
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FIG. 5: The calculated spatial temperature distribution of a 
para junction ([1,4]BDT) with Ti = 3257f and T 2 = 275K, 
measured by an SThM scanned 3.5A above the plane of the 
carbon nuclei. The Pt SThM tip is assumed to be atomically 
sharp, with an effective blackbody surface area equal to that 
of a 50nm sphere with an emissivity of 0.1. The sulfur linker 
atoms are indicated schematically by red and blue circles in- 
dicating the contacts to the hot and cold leads, respectively. 
Quantum oscillations of the temperature are clearly visible in 
the vicinity of the molecule, which can be explained in terms 
of the Kekule contributing structures shown in Fig. [7] 



ent from that of its nearest neighbors: orbitals 2 and 6 
are hot, orbitals 3 and 5 are cold, while orbitals 1 and 
4, directly connected to the hot and cold electrodes, re- 
spectively, have intermediate temperatures. For weaker 
thermal coupling of the probe to the ambient environ- 
ment, thermal oscillations are observed at even larger 
length scales, while stronger coupling to the environment 
reduces resolution of these quantum effects. 

Quantum oscillations of the temperature in a nanos- 
tructurc subject to a temperature gradient are a ther- 
mal analogue of the voltage oscillations predicted by 
Buttikcr 13 in a quantum system with an electrical 
bias. Similar temperature oscillations in one-dimensional 
quantum systems were first predicted by Dubi and Di 
Ventral 

In order to understand the temperature oscillations 
shown in Fig. [5j it is useful to consider the thermal con- 
ductance hp\ between the probe and the hot electrode. 
Fig. [6]^ a) indicates that k p p is large when the probe is 
in the ortho or para configuration relative to electrode 
/?, as well as when it is proximal to the 7r-orbital di- 
rectly coupled to the electrode. However, k p p is nearly 
zero when the probe is in the mcta configuration rela- 
tive to electrode /?. The three-terminal correction to the 
thermal conductance Ak p i is plotted in Fig. Efb) , which 
indicates that three-terminal thermoelectric effects lead 
to a sizable relative correction to the thermal conduc- 
tance between the probe and the electrode in the meta 
configuration, but are otherwise small. 

Eq. (JSj> for the local temperature can be rewritten in 
the following instructive form: 



T P — Tq 



(k p i - K p2 )AT 

k p l + K p 2 + KpQ ' 



(14) 



where Ti = T + AT/2 and T 2 =T - AT/2. Thus, when 
Kpi 3> k P 2, the probe will measure a temperature near 
T\ . and vice versa when k p i <C k P 2, provided the coupling 
to the environment is not too large. On the other hand, 
the probe will measure a temperature near To if h p \ ~ 
k P 2- Comparing Figs. [5] and El one sees that 7r-orbitals 2 
and 6 are hot because when the SThM is coupled locally 
to them, it is in an ortho configuration relative to the 
hot electrode and a meta configuration relative to the 
cold electrode, and Kortho 3> £ mcta . Orbitals 3 and 5 are 
cold by symmetry. On the other hand, orbitals 1 and 4 
have intermediate temperatures since k p \ ~ h P 2 when the 
SThM is in a para configuration relative to one electrode 
and proximal to the other. 

The quantum oscillations in the temperature shown 
in Fig. [5] can also be understood as consequences of the 
rules of covalence in conjugated systems. Fig.[7]shows the 
Kekule contributing structures illustrating charge trans- 
fer from an electrode E to a benzene molecule. Consid- 
ering both hot and cold electrodes, the rules of covalence 
dictate that electrons from the hot electrode are available 
to tunnel onto the temperature probe when it is coupled 
locally to orbitals 2, 4, or 6, while electrons from the cold 
electrode are available to tunnel when the probe is near 
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FIG. 6: (a) Thermal conductance k p \ between lead 1 (hot) and the temperature probe p. k p \ is largest when p is in an ortho 
or para configuration relative to the hot electrode, or proximal to it (near orbitals 1, 2, 4, or 6), and smallest when it is in a 
meta configuration relative to the hot electrode (near orbitals 3 and 5). (b) The difference between the three- and two-terminal 
thermal conductances k p i — k p i, showing the largest errors (1.03 x 10 _6 k;o) occur where k p i is small (near sites 3 and 5), 
indicating large relative errors if a two-terminal formulation were used. Here Ko = (n 2 /3)k%T /h = 2.835 x W~ 4 eV/sK. 




FIG. 7: Kekule contributing structures illustrating charge 
transfer from an electrode E onto benzene. The rules of co- 
valence in conjugated systems dictate that electrons from an 
electrode are available to tunnel onto the temperature probe 
when it is coupled locally to the molecule in an ortho or para 
configuration relative to the electrode. 
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FIG. 8: The calculated spatial temperature distribution of an 
ortho junction ([1,2]BDT) measured under the same condi- 
tions described in Fig. [5] The quantum temperature oscilla- 
tions can be explained in terms of the Kekule contributing 
structures shown in Fig. [7] 



orbitals 1, 3, or 5. In addition, electrons from the hot 
electrode are available to tunnel onto the temperature 
probe when it is near orbital 1, which is proximal to the 
hot electrode, while electrons from the cold electrode are 
available to tunnel onto the temperature probe when it 
is near orbital 4, which is proximal to the cold electrode. 
Orbitals 2 and 6 thus appear hot, orbitals 3 and 5 appear 
cold, while orbitals 1 and 4 should exhibit intermediate 
temperatures by this argument. 

The calculated temperature distribution of an ortho 
junction is shown in Fig. |8l measured under the same 
conditions discussed in Fig. [5J The Kekule contribut- 
ing structures illustrated in Fig. [7] dictate that lone pairs 
from the hot electrode may tunnel to orbitals 2, 4, or 



6, and lone pairs from the cold electrode may tunnel to 
orbitals 1, 3, or 5. Taking into account that orbitals 1 
and 2 are also proximal to the linker groups binding the 
molecule to the hot and cold electrodes, respectively, and 
thus exhibit intermediate temperatures, the rules of co- 
valence dictate that orbitals 4 and 6 appear hot, while 
orbitals 3 and 5 appear cold, in complete agreement with 
the calculated temperature distribution. 

For the para and ortho junctions shown in Figs. [5] and 
[HI the rules of covalence act essentially like a Maxwell 
demon, in that they selectively permit electrons from the 
hot or cold reservoir to tunnel onto the probe when it is 
at specific locations near the molecule, and block elec- 
trons from the other reservoir. The question might arise 
whether the actions of this Maxwell demon could lead 
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FIG. 9: The calculated temperature distribution of a meta 
junction ([1,3]BDT) measured under the same conditions de- 
scribed in Fig. [5] In this case, the temperature distribution 
exhibits well-defined "bond temperatures" arising from off- 
diagonal contributions to the thermal transport, in addition 
to the "orbital temperatures" of orbitals 1 and 3 proximal to 
the hot and cold electrodes, respectively. 



to a violation of the second law of thermodynamics, as 
Maxwell originally hypothesized. However, in this case 
there is no violation of the second law, because electrons 
within the molecule "remember" which electrode they 
came from. There is no "mixing" of the hot and cold 
electrons in the absence of inelastic scattering, which 
is strongly suppressed compared to elastic processes in 
these junctions at room temperature. 

The calculated temperature distribution of a meta 
BDT junction is shown in Fig. [9j measured under the 
same conditions as in Figs. [5] and [H In this case, the 
temperature distribution is more complicated, exhibiting 
both well-defined "orbital temperatures" (1 and 3) and 
"bond temperatures" (4-5 and 5-6 bonds). The tem- 
peratures of orbitals 1 and 3 can be explained by the 
arguments given above, while the rules of covalence illus- 
trated in Fig. [7] indicate that neither electrons from the 
hot electrode nor the cold electrode can reach orbital 5, 
so that its temperature is indeterminate. Near orbital 5, 
off-diagonal contributions to the transmission dominate 
due to the suppression of transmission in the meta con- 
figuration. The para transmission amplitude interferes 
constructively with the small but nonzero meta transmis- 
sion amplitude, while the ortho transmission amplitude 
interferes destructively with the meta transmission am- 
plitude, so that the 4-5 bond appears hot while the 5-6 
bond appears cold. 



B. [l,6]Hexatrienedithiol junction 

The calculated temperature distribution of a [l,6]hexa- 
trienedithiol junction composed of a thioloated 6-site lin- 
ear molecule (hexatriene) covalently bonded to two gold 
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FIG. 10: The calculated temperature distribution for a Au- 
[l,6]hexatrienedithiol-Au molecular junction, measured under 
the same conditions discussed in Fig. [5] The observed temper- 
ature oscillations can be explained using the resonance con- 
tributing structures shown in Fig.rn] 




FIG. 11: Resonance contributing structures for a hexatriene 
junction. For a given electrode E, the lone pair can only 
occupy every other 7r-orbital, giving rise to the alternating 
hot, cold temperature profile shown in Fig. 1101 

electrodes is shown in Fig. [TO] The conditions of the 
temperature measurement are the same as described in 
Fig. [5] Quasi-one-dimensional temperature oscillations 
are clearly observable along the length of the molecular 
wire, consistent with the prediction of Ref. OH The reso- 
nance contributing structures describing electron transfer 
from an electrode E onto the molecule are shown in Fig. 
ITT1 As in the case of the para and ortho-BDT junctions, 
the rules of covalence are unambiguous, and predict alter- 
nating hot and cold temperatures for the 7r-orbitals along 
the length of the molecule, with intermediate tempera- 
tures for the end orbitals proximal to the two electrodes, 
consistent with the calculated temperature distribution. 

C. [2,7]Pyrenedithiol junction 

As a final example, we consider the effects of finite spa- 
tial resolution on the measured temperature distribution 
of a polycyclic [2,7]pyrenedithiol-Au junction. The tem- 
perature distribution was calculated for three different 
values of the SThM spatial resolution in the three pan- 
els of Fig. [T^l The leftmost panel shows the maximum 
spatial resolution under the specifications of a hypothcti- 
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FIG. 12: The calculated spatial temperature distribution for a Au-[2,7]pyrenedithiol-Au junction with T± — 325K and T2 = 
275K, measured at a height of 3.5A above the plane of the molecule, for three different values of SThM spatial resolution. 
The leftmost panel shows the maximum spatial resolution under the specifications of a hypothetical SThM described in Fig. 
[5] The middle and rightmost panels show the measured temperature distributions with reduced spatial resolution obtained 
by convolving the tip-sample tunnel-coupling T p with Gaussian distributions with standard deviations of a = 0.5A and 1.0A, 
respectively (full resolution corresponds to a = 0). In the rightmost panel, quantum temperature oscillations are no longer 
resolved, and the temperature distribution resembles a thermal dipole. Similar results were obtained for a variety of molecular 
junctions, suggesting that a classical temperature distribution consistent with Fourier's law of heat conduction emerges due to 
measurements with limited spatial resolution. 



cal SThM described in Fig. [5l The middle and rightmost 
panels show the measured temperature distributions with 
reduced spatial resolution obtained by convolving the tip- 
sample tunnel-coupling T p with Gaussian distributions 
with standard deviations of a = 0.5A and I.OA, respec- 
tively (full resolution corresponds to a = 0). Many-body 
transport calculations for this larger molecule are cur- 
rently computationally intractable, so we have utilized 
Hiickel theory to describe the molecular electronic struc- 
ture, as discussed in the Supporting Information. 

The left-most panel of Fig. [12] shows a complex inter- 
ference pattern of hot and cold regions with a symmetry 
that mimics the junction itself (in this case with two mir- 
ror axes). More complex molecules such as this highlight 
the "proximity effect" whereby the flow of heat from a 
given electrode to the orbitals in its vicinity is enhanced, 
so the molecule is generally warmer near the hot lead and 
cooler near the cold lead. We mention that, although 
more tedious, the Kekule contributing structures^ can 
be used to understand the pattern of temperature varia- 
tions in this molecule as well. 

Focusing on the middle and right-most panels of 
Fig. 1121 we find an immediate consequence of the proxim- 
ity effect as a is increased: non-monotonic temperature 
variations due to quantum interference are washed out, 
and the underlying temperature gradient appears. In 
the rightmost panel, quantum temperature oscillations 
are no longer resolved, and the temperature distribu- 
tion resembles a thermal dipole. Similar results were 
obtained for a variety of molecular junctions, suggest- 
ing that a classical temperature distribution consistent 
with Fourier's law of heat conduction emerges when the 



temperature is measured with limited spatial resolution. 

The transition from microscopic quantum tempera- 
ture oscillations to macroscopic diffusive behavior and 
Fourier's law is still poorly understood^ It has been ar- 
gued that Fourier's law is recovered in systems with suf- 
ficient dephasing^ or disorder £ However, we have seen 
that in conjugated organic molecules, the quantum tem- 
perature oscillations arc intimately connected to the rules 
of covalence describing the 7r-bonds of the molecule. De- 
phasing (or disorder) sufficient to wash out the tempera- 
ture oscillations would thus necessarilly sever the 7r-bonds 
and dissociate the molecule. Since the molecules studied 
in this article arc stable at room temperature, we know 
that such strong dephasing cannot be present. Thus we 
predict that quantum temperature oscillations will be ob- 
served in molecular junctions if temperature measure- 
ments with sufficient spatial resolution are performed, 
and that Fourier's law is a consequence of coarse-graining 
due to finite spatial resolution. 



IV. CONCLUSIONS 

We have proposed a physically motivated and mathe- 
matically rigorous definition of an electron thermometer 
as an electron reservoir coupled locally to and in both 
thermal and electrical equilibrium with the system be- 
ing measured [cf. Eq. (JTJ)]. This definition is valid under 
general nonequilibrium conditions with arbitrary thermal 
and/or electric bias. Based on this definition, we have de- 
veloped a realistic model of an atomic-resolution SThM 
operating in the tunneling regime in ultrahigh vacuum, 
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where the resolution of temperature measurements is lim- 
ited by the radiative thermal coupling of the probe to the 
ambient black-body environment. 

We used this model of an atomic-resolution SThM to 
investigate the nonequilibrium temperature distributions 
of a variety of single-molecule junctions subject to ther- 
mal gradients. Quantum oscillations of the local tem- 
perature that can be observed using a SThM with suf- 
ficiently high resolution are predicted. We show that in 
many cases, these quantum temperature oscillations may 
be understood straightforwardly in terms of the rules of 
covalence describing bonding in 7r-electron systems. As 
such, these oscillations are predicted to be extremely ro- 
bust, insensitive to dephasing or disorder that is insuffi- 
cient to dissociate the molecule. Instead, we show that 
such quantum interference effects arc washed out if the 
spatial resolution of the SThM is insufficient to observe 
them, and that the temperature distribution then ap- 
proaches that expected based on Fourier's classical law 
of heat conduction. Message: The temperature oscilla- 
tions are really there, if you look closely enough! 

One may wonder whether it is meaningful to define a 
temperature that varies significantly from place to place 
at the atomic scale. Since temperature is related to mean 
thermal energy, does a variation of temperature on a scale 
comparable to the de Broglie wavelength not violate the 
uncertainty principle? Our answer to such questions is a 



pragmatic one: By definition, temperature is that which 
is measured by a thermometer, and the position of a 
thermometer can certainly be controlled with subatomic 
precision using standard scanning probe techniques. We 
should also emphasize that our proposed thermometer 
measures the electron temperature, which may be largely 
decoupled from the lattice temperature in nanoscale junc- 
tions. 

Finally, let us return to the theme of the title of this 
article. We have shown that in a molecular junction con- 
taining a conjugated organic molecule, the rules of cova- 
lence act essentially like a Maxwell demon, in that they 
selectively permit electrons from the hot or cold reservoir 
to tunnel onto the probe when it is at specific locations 
near the molecule, and block electrons from the other 
reservoir. The question might arise whether the actions 
of this Maxwell demon could lead to a violation of the 
second law of thermodynamics, as Maxwell originally hy- 
pothesized. However, in this case there is no violation of 
the second law, because electrons within the molecule 
"remember" which electrode they came from. There is 
no "mixing" of the hot and cold electrons in the absence 
of inelastic scattering. And we have argued that dephas- 
ing due to inelastic scattering is insufficient to perturb 
this particular embodiment of Maxwell's demon without 
dissociating the molecule itself. 
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